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The proton temperature anisotropy in the solar wind is known to be constrained by the theo- 
retical thresholds for pressure anisotropy-driven instabilities. Here we use approximately 1 million 
independent measurements of gyroscale magnetic fluctuations in the solar wind to show for the 
first time that these fluctuations are enhanced along the temperature anisotropy thresholds of the 
mirror, proton oblique firehose, and ion cyclotron instabilities. In addition, the measured magnetic 
compressibility is enhanced at high plasma beta (f3\\ > 1) along the mirror instability threshold 
but small elsewhere, consistent with expectations of the mirror mode. The power in this frequency 
(the 'dissipation') range is often considered to be driven by the solar wind turbulent cascade, an 
interpretation which should be qualified in light of the present results. In particular, we show that 
the short wavelength magnetic fluctuation power is a strong function of collisionality, which relaxes 
the temperature anisotropy away from the instability conditions and reduces correspondingly the 
fluctuation power. 
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The physical processes that regulate the expansion of 
the super-Alfvenic solar wind into space include adia- 
batic particle motion, plasma instabilities, and binary 
particle collisions. As the wind expands, plasma den- 
sity n p and magnetic field \B\ decrease radially. The 
Chew-Goldberger-Low (CGL) relations [1 predict that 
the plasma ions should become anisotropic in the sense 
of Tj| > T± if the particle motion is adiabatic and col- 
lisionless; here T is the ion temperature parallel and 
perpendicular to the background magnetic field. How- 
ever, Coulomb collisions and pressure-anisotropy insta- 
bilities act to pitch-angle scatter the plasma back towards 
isotropy [2]. At 1 AU, the most probable value of the pro- 
ton temperature anisotropy is T±/T\\ « 0.89 (see Figure 
1, top panel below). If CGL were valid, this would imply 
a proton temperature anisotropy of T±/T\\ > 200 at 5 
solar radii. 

The same pressure-anisotropy instabilities that oper- 
ate in the solar wind are believed to operate in other 
low-collisionality astrophysical plasmas, including clus- 
ters of galaxies [3] and some accretion disks onto black 
holes |H [3]. In the latter environment, these instabilities 
not only modify the thermodynamics of the plasma (as in 
the solar wind), but they also play a crucial dynamical 
role, regulating the anisotropic stress that helps trans- 
port angular momentum, allowing accretion to proceed. 

Growth of ion temperature anisotropy instabilities has 
been studied in a (relatively collisional) laboratory de- 



* Physics Department and Space Sciences Laboratory, University of 
California, Berkeley 

t Harvard-Smithsonian Center for Astrophysics, Cambridge 
■'•Department of Physics and Astronomy, University of Iowa, Iowa 
City 

§ Physics Department and Astronomy Department, University of 
California, Berkeley 

^ Space Sciences Laboratory, University of California, Berkeley 



vice, where isotropization and magnetic fluctuations were 
observed corresponding to the Alfven ion cyclotron in- 
stability [6 . Similar results were obtained in the solar 
wind at 1 AU [7] (for T±/T\\ > and solar wind speeds 
greater than 600 km/s), suggesting that the proton cy- 
clotron instability plays a role. Both mirror mode and 
ion cyclotron anisotropy instabilities appear to be active 
in the terrestrial magnetosheath, inferred from observed 
constraints on the temperature anisotropy [8] [9] . 

More recently, Kasper et al. [10] used measurements 
by the SWE instrument [TT] on the Wind spacecraft to 
show that the solar wind proton temperature anisotropy 
T±/Tn is constrained with respect to proton parallel beta 
(3\\(= n p ki)T\\ / B 2 ) in a way that is consistent with expec- 
tations of the proton firehose instability (for T_\_/T\\ < 1). 
Hellinger et al. (2006) compared linear instability calcu- 
lations with the measured anisotropics (from this same 
dataset) and found that the constraints on the observed 
anisotropics were best described by the mirror instabil- 
ity (in contrast to previous results [7]) for the case of 
perpendicular anisotropy (Tj_/Tj| > 1) and the oblique 
firehose instability for parallel anisotropy (Tj_/Tj| < 1). 
An instability threshold can be derived by calculating 
marginal stability (in practice 7 = 10 -3 uo C i) for values of 
[12]. 

T± _ a 

+ (/?||-/3 ) b { ' 

Equation (1) generalizes previous results [71 [13] and 
constrains the Wind/SWE observations for values of a, 6, 
and Pq corresponding to the mirror and oblique firehose 
instabilities (coefficients are given below). 

Here we use experiments on the Wind spacecraft to 
show the first direct correspondence of the measured 
wave power to the anisotropy-driven instabilities. Wind 
was launched in 1994 and has spent long intervals in the 
solar wind, often near the first Lagrange point upstream 
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of the Earth. The data set spans the time interval Nov. 
21, 1994 to Nov. 12, 2004 and includes 1,026,112 inde- 
pendent measurements of solar wind plasma and mag- 
netic field. Proton density n p , velocity v SWl and tem- 
perature are measured by the Faraday cup instrument 
of SWE (Solar Wind Experiment) [IT]. Both the par- 
allel proton temperature Tji and perpendicular temper- 
ature T± are computed by comparison with the average 
magnetic field direction. The Magnetic Field Investiga- 
tion (MFI) is used to measure the solar wind magnetic 
field [14 at either 22 vectors/s or 11 vectors/s depend- 
ing on telemetry mode and averaged to 3 second inter- 
vals; we denote this 3 second average data as B. The 
vector rms fluctuation field during the 3 second average 
interval is denoted as SB. Each fluctuation field mea- 
surement SB is rotated into a coordinate system defined 
by the average field direction £?, so that we have both 
the compressive component 5B\\ and the shear compo- 
nent 5B ± (= (SB 2 ±1 + 5B\ a fl 2 ). We then define the 
magnetic compressibility as SB 2 /(SB 2 + 5B\) [15 j. The 
magnetic field is an rms measurement over the bandwidth 
Af = (0.3 - 5.5) Hz, or (0.3 - 11) Hz in high TM mode. 
Since the magnetic fluctuation spectrum falls as f~ 5 ^ 3 (or 
steeper) at these frequencies, the power in this bandwidth 
is dominated by the lowest measured frequency fo ~ 0.3 
Hz. Since it is believed that k± ^> ku [16], the natu- 
ral frequencies of the turbulence are all much lower (i.e. 
fo ^> fd) an d Taylor's hypothesis applies (uj = kv sw ). 
Therefore the measured power at /o ~ 0.3 Hz corresponds 
to power at wavenumber kp { w (f / ' f C i){v th /v sw ). The 
distribution (in our data) of this parameter is sharply 
peaked at kpi « 0.56 with a half- width of 0.32; therefore 
these measurements correspond to magnetic fluctuation 
power at kpi « 0.56 ± 0.32. 

The upper panel of Figure 1 shows the distribution of 
measurements of proton temperature anisotropy (T±/T\\) 
against parallel proton plasma beta f3\\] this distribu- 
tion shows the striking signatures of the regulation of 
the anisotropy by instabilities. These proton measure- 
ments are a subset of those used by Kasper et al.[10 ; 
here we include only proton measurements at times when 
there also exist good measurements of SB. Dotted lines 
(on all four panels) show the instability thresholds for 
the mirror instability at T±/T\\ > 1 (equation 1 with 
(a,b,/3 ) = (0.77,0.76,-0.016)) and the oblique firehose 
instability at T_\_/T\\ < l(equation 1 with (a, /?o) = 
(-1.4,1.0,-0.11)) [12]. The second panel of Figure 1 
shows the average measured amplitude of the magnetic 
fluctuations in the space of (/3|i,Tj_/Tji), as in 

the upper panel. A general trend of enhanced fluctua- 
tions with larger (3\\ is clearly visible. Furthermore, the 
fluctuation amplitude is enhanced along both the mir- 
ror /IC (for T±/Tn > 1) and oblique firehose (T±/T\\ < 1) 
thresholds. This is strong evidence that these instabil- 
ities are excited here and act to isotropize the plasma; 
this is the principal result of this study. The third 
panel of Figure 1 shows the magnetic compressibility 
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FIG. 1: The distribution of proton temperature anisotropy 
(Xj_/Tj|) measurements with respect to the parallel plasma 
beta (3\\ (upper panel) is constrained by the oblique proton 
firehose instability threshold (lower dotted line) and the mir- 
ror instability threshold (upper dotted line). In the second 
panel, the magnetic fluctuation amplitude |&B|/|.B| is shown 
to be enhanced along the instability thresholds and overall at 
high (3\\ where the thresholds converge. The third panel shows 
magnetic compressibility SB^/(SB^-\-SB^), which is enhanced 
at high f3\\ (>1) along the mirror threshold, as expected for the 
mirror instability. The lower panel shows the 'collisional age', 
which is largest around T±/T\\ « 1 suggesting that isotropy 
results largely from Coulomb collisions. Anisotropic plasma 
is relatively collisionless. 
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5B 2 /(SB 2 + 5B±), which is enhanced to values of ~ 0.3 
for compressive solar wind (Tj_/Tj| > 1) with f3\\ > 1, as 
would be expected for the mirror instability. The com- 
pressibility becomes smaller for f3\\ < 1, which is con- 
sistent with the Alfven ion cyclotron mode, however the 
power continues to be bounded by the mirror threshold. 
Linear mirror instability calculations [15] for T±/T\\ > 1 
predict values of the magnetic compressibility between 
0.8 and 1, therefore our measurements suggest a mixture 
of waves. Furthermore, it is interesting to note that the 
typical value of the magnetic compressibility away from 
the thresholds is small ~ 0.1, in contrast to theories which 
suggest evolution of the MHD cascade to the (compres- 
sive) magnetosonic/whistler branch at short wavelengths 
[T7] . The bottom panel of Figure 1 shows the average 
collisional 'age' in each (f3^T±/Tu) bin; the collisional 
age r is defined as t — i^pp 

L/v sw , the Coulomb proton- 
proton collision frequency v pp multiplied by the transit 
time from the Sun to 1 AU and is an estimate of the 
number of binary collisions in each plasma parcel during 
transit from the Sun to the spacecraft. It is interesting, 
however obvious, that the more collisional plasma is more 
isotropic; away from T±/T\\ « 1 the plasma is relatively 
collisionless. It has been shown recently [18] that colli- 
sional age organizes solar wind instabilities better than 
the traditional distinction of 'fast' and 'slow' wind. 
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FIG. 2: The magnitude of magnetic fluctuations \5B\ aver- 
aged into bins of T±/T\\ vs f3\\. Enhanced power exists well 
away from the thresholds, as expected. The regions of en- 
hanced SB corresponds to the enhanced proton heating in 
Liu et al (2006). 

Figure 2 shows the magnetic fluctuations data \5B\ un- 
normalized by Bo. Linear instability thresholds associate 
with the mirror, ion cyclotron (AIC) , and oblique fire- 
hose instabilities [12] are overlaid. It is interesting to note 
that, as found by Hellinger et al [12], the oblique firehose 
and the mirror instabilities appear to constrain the ob- 
served distribution of data, not the ion cyclotron nor par- 
allel firehose instabilities (in spite of a larger growthrate); 



this may be because both the mirror and oblique fire- 
hose are non-propagating instabilities. The regions of 
enhanced magnetic fluctuations, near the mirror and 
oblique firehose thresholds also correspond to measure- 
ments of enhanced proton temperature published else- 
where [19] , It is unclear if this indicates plasma heating 
due to anisotropy instabilities, in addition to pitch-angle 
scattering, or if the 'younger' (less collisional) plasma is 
merely hotter than average. 
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FIG. 3: Magnetic fluctuation amplitude \SB\ 2 as a function 
of the collisional age; the white dots are the most likely value 
of 1 5B | 2 in each age bin. Magnetic fluctuations near the pro- 
ton gyroradius (kp « 1/2) are suppressed in more collisional 
plasma. Coulomb collisions maintain the isotropy of the pro- 
tons, which then remain far from the instability thresholds. 
Note that this corresponds to a factor of 100 suppression of 
magnetic power 5B 2 over the full range of collisionality. 

Figure 3 shows histograms of the fluctuation ampli- 
tude squared \SB\ 2 in bins of collisional age; the white 
dots show the most probable value. The overall mag- 
netic fluctuation power SB 2 is a function of the colli- 
sional age, with the magnetic power weaker by a fac- 
tor of ~100 for more collisional plasma. This effect is a 
proxy for the temperature anisotropy: collisional plasma 
is more isotropic and, therefore, further from the instabil- 
ity thresholds. This underscores the important point that 
the power spectral density (PSD) of magnetic fluctuation 
power near 1 Hz in the solar wind is modified by these 
local instabilities. Previous studies of the PSD of short 
wavelength solar wind turbulence have not accounted for 
this and should be reexamined [2Q1 EE] E2 E3] [24]. If 
the | SB | 2 values in Figure 3 are divided by the measure- 
ment bandwidth (approximately 5-10 Hz depending on 
sample rate), they can be compared to power spectral 
density (PSD) measurements published previously (over 
the bandwidth of 0.3 to 11 Hz), noting that the power 
is dominated by the amplitude at the lowest frequencies 
(0.3 Hz). A log- log fit to the most probable values in 
Figure 3 gives a relationship \SB\ 2 w 10 _4 r -0,57 , where 
T — v vv L/v sw is the collisional age. 
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FIG. 4: Wavelet power spectra of magnetic fluctuations at 
three different time intervals corresponding to perpendicu- 
lar anisotropy T±/T\\ = 2.2, /3\\ = 0.2 (red trace), parallel 
anisotropy T±/T\\ = 0.5, /3\\ = 1.9 (blue trace) and isotropic 
protons T±/T\\ = 1, /3\\ = 0.7 (black trace). The spectra are 
plotted against kp and have been scaled to the same average 
value over the interval kp G (0.2, 1.0), as described in the text. 
A solid black line with scaling k~ 5 ^ 3 is shown. The interval 
above kp « 1/2 corresponds to the statistical data shown in 
Figures 1 and 2. 



Figure 4 shows wavelet power spectral density for 3 
distinct time intervals; these spectra are shown plotted 
against kpi assuming Taylor's hypothesis (lj = kv sw ) 
and using the local values of v sw ,vth, and uo C i and scal- 
ing the spectra to have similar power over the range 
kpi G (0.2,1.0). The black trace in Figure 4 shows 
the magnetic power from an interval of isotropic plasma 
(Xj_/Tj|, /3||) ~ (1,0.7), the red trace is data from an inter- 
val of perpendicular anisotropy (Tj_/Tj| , (3\\) « (2.2,0.2), 
and the blue trace from an interval of parallel anisotropy 
(Tj_/Tj|,/3||) w (0.5,1.9). Although wavenumbers below 
kpi w 0.2 are aliased by the Morlet wavelet transform 
and are unreliable. The inert ial range, where the power 
is proportional to k~ 5 ^ 3 (shown by a solid black line in 
Figure 4), can be made to agree, the nature of the short 
wavelength spectrum above kpi « 1 varies between spec- 
tra. While this result is not statistical, nor is it surprising 
given the clear organization of the magnetic fluctuation 
power in Figure 1, it begs a reexamination of 'dissipation 
range' turbulence. Future observational studies of tur- 
bulent dissipation and secondary (kinetic Alfven wave) 
cascades should attempt to understand the role of these 
local instabilities, which may confuse the interpretation 
of solar wind magnetic power spectra. 
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